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ABSTRACT 

The effect of surface roughness on the brightness temperature of a moist 
terrain has been studied through the modification of Fresnel reflection co^i- 
cient and using the radiative transfer equation. Model calculations are in good 
qualitative agreement with the observed dependence of the brightness tempera* 
ture on the moisture content in the surface layer. 

Introduction 

There have been several recent attempts to calculate the brightness temper- 
ature using the radiative transfer model in solid medium (e. g. , Stogryn, 1970; 

England, 1974; Tsang and Kong, 1975). These attempts were primarily directed 
towards the understanding of microwave emission from snow fields. Wilheit 
(1975) and Burke and Paris (1975) have performed model calculations for the 
microwave emission from soil assuming the medium to be a layered structure. 

One common assumption in all these calculations is that the layer separating 
the medium emitting the radiation is a geometrically smooth surface. For micro- 
wave emission from snow fie .ds or deserts, this assumption may have some justi- 
fication but for a natural terrain or an agricultural field, this assumption is highly 



(piestionable. In fact, recent observation by Schmugge et. al. (1976a) of the 
variation of the brightness temi^rature with the soil moisture in the surface 

layer show deviation from the radiative transfer model calculations. The 
purpose of this paper is to show the ramification of surface roughness on the 
brightness temperature. The smooth surface Fresnel reflection co^ficient 
has been modified by assuming a Gaussian height distribution such that at all 
points on the surface the slope is less than unity and the radii of curvature 
are greater than the radiation wavelength (Barrick 1970). The present study 
differs from the one i^rformed by Wu and Fung (1972) for the microwave emis- 
sion from sea. They used a composite-surface scattering theory to include the 
effect of surface roughness but did not use the radiative transfer equation to 
calculate the apparent temperature. Although the composite-surface model 
seems more appropriate for a rough surface, it is computationally prohibitive 
without further approximations and also introduces several parameters to 
characterize the surface. The model used in this paper is simpler in that it 
has only one parameter to characterize the surface, namely, the standard 
deviation in height distribution. It is shown that the present model can provide 
a qualitative account of the dependence of the brightness temperature on the 
surface layer soil moisture. 

Formulation 

To describe the microwave emission from the soil, we will consider the 
radiative transfer equation (Chandrasekhar, 1960): 
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where I is the Intensity propagating in the direction u, K. is the t(^ extinction 
per unit lei^^ and S is the source term describing the contribution to the in- 
tensity due to scattering and due to the continuum thermal emission of the soil. 

In principle one should study this equation in conjuiu;tion with the equation de- 
scribing the heating of the soil. It is this latter equation which will provide in- 
formation ^out the thermal part of the source term. In this paper we will decouple 
these two ^nations in the sense that we will consider a given temperature 
distrilHition. 

To solve the radiative transfer equation we will ccmsider a semi-infinite 
medium with dej^h dependent temperature and moisture distrilxitions. Since 
soil is a highly absorbing material (i.e., large imaginary part of dielectric con- 
stant), to a good apprmcimation, the brightness temperature Tg or the tempera- 
ture equivalent of the intensity emerging from the soil will be ctetermined by its 
internal temperature distribution T(z). By integrating eqn (1) with the source 
term as the temperature distribution of the soil one can write (Burl^ and Paris 
(1975)) 

dz'/cos i>^dz 

where rp is the soil surface reflectivity for polarization p. The absorption length 
and the angle cos 6 are related to the angle of observation 6 ^ and the proper- 
ties of the medium as follows (Born and Wolf, 1975). 

If the complex refractive index of the soil is written as: 


(1 - 


r ) 


r 


T(z) (Kg/cos 9 ) exp 



N = n(l + iK)= /T 
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then 


cos 6 = 


n q (cos y - K sin y) 

[sin^ dQ + q* (cos 7 - k sin 


- — nq (^< cos y + sin y) 

X. 


where q and 7 are defined by 


q2 cos 27=1 ^ sin^ 0 , 


2^2 


n ^ (1 + 


q2 sin 27= ^ ^ sin^ 0 ^ 

n^ (1 + 

and is the free space wavelength of the radiation and e is the complex dielectric 
constant of the medium which depends iq)on the wavelength and the moisture of the 
soil. Thus from the knowledge of the temperature and the moisture distributions, 
one can calculate the brightness temperature by performing the integral in eqp 
(2). This approach to the calculation of the brightness temperature is similar to 
the models developed by Wilheit (1975) and Burke and Paris (1975). The models 
txeat the medium as consisting of layers with constant moisture and temperature. 
We will now consider the effect of the soil surface roughness on the br ipht- 
nesB temperature. A first principle account of this effect is extremely difficult 
and far from being of practical use. In this paper we will consider the following 
simpler approach which requires modification of eqn (2). 

For a vertically stratified medium, the electric field within the medium can 
be represented by 

E = Eq exp [i (k^ X + 2 + k^ y)J (3) 
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k . k , and k, are the x, y and z components of the complex wave vector 

A y z 

with z-axis being pei^ndicular to soil surface. This wave will now be con- 
sidered to be scattered by the Interface of the medium and the free space. 


It has been shown by Tolstoy and Clay (1966) that if this interface is a sta- 
tistically rough surface such that there is no correlation between the amplitudes 
of the waves scattered by two points on the surface, then the scatter^ intrasity 
can be obtained by the absolute square of the average scattered amplitude. It has 
further been shown that if represents the scattered amplitude by a perfectly 
smooth and perfectly reflecting surface, then the average amplitude that will be 
scattered by a rough surface is given by 


<E.o> = *0 f 


ik 


dz 


(4) 


where W(z) is the height distribution of the surface and Ro is the reflection 

coefficient of a smooth surface. The result for a smooth surface can be 

detained trivially by identifying the distribution with a delta function. A typical 
rough surface corresponds to identifying the spectrum with a Gaussian distribu-f 

tion of zero mean and variance a - : 


\V(z) = — L— .exp [-zV2cr^] 

a V " 2 77 


(5) 


For this spectrum, the average amplitude is given by 

<E.o> (6) 

By demanding the continuity of phase at the interface we can write for a 



wave in x-z plane 
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Since within the medium we have: 




so that 




e - sin-^ 


Thus the scattered intensity is (Stained from equation (6): 
I, = ij IRq exp [-h (Re t - Sin^ 6*^)] 


(7) 


(8) 


where Re t represents the real part of dielectric constant and the roughness 
parameter h is given by 




(9) 


From this result, one can stipulate that the gross effect of the surface rough' 
ness on the scattered intensity can be incorporated by modifying the smooth sur- 
face reflectivity r„p (= |Ro | ^ ) as 

•■p = '■qp ^ (10) 

whore the subscript p designates the polarization. 

Note that this result differs from the result obtained in the radar cross- 
section analysis (Barrick, 1970): 

'■p - •'op ‘‘’‘P (- *o) (11) 

because in the present case the energy is incident on the interface from within 
the soil as opposed to the radar case where the energy is incident from air. 
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Results and PiscusBion 

Experimental results reported in this paper were obtained chiring aircraft 
flights over the Phoenix, Arizona area and the Imperial Valley of California dur- 
ing March 1972 and February 1973 (Schmugge et. al. 1976a) and flights over only 
the Phoenix area during March 1975 (Schmugge 1976b). The surface roughness 
characteristics were those resulting from the agricultural practices of the two 
areas. The dominant method of irrigation is the flooded furrow. The furrow 
separation was about one meter aiul the furrow hei^t was about 20 cm. Super- 
imposed on these corrugations were clods, which were gemrally less than 5 cm. 

Calculations were performed ^ wavelengths of 1. 55, and 21 cm using 
moisture and temperature profiles observed at the U. S. Water Conservaticm 
Laboratory at Phoenix in March 1971 (see Schmugge et. al. 1976a). These 
profiles can be assumed to be reasonable estimates for the actual situation for 
the data obtained in March 1972. For the data obtained in February 1973, the 
surface temperature was found to be about 15K lower than that observed during 
March. The February temperature profiles were obtained from the observed 
March temperature profiles by adjusting the gradient to fit the deserved surface 
temperature for February data and assuming two profiles to be equal at about 
30 cm. For the data obtained in March 1975, the calculations were performed 
usii^ the actual ground measurement of the temperature and moisture profiles 
at each site (Birchard 1975). The dependence of the real and the imaginary part 
of the dielectric constant on the weight percent of soil moisture for different 
wavelengths used in this calculation were obtained by Schmugge et. al. (1976a) 
through the linear regression fit of the data. The results for 1.55 and 21 cm 
wavelength is shown in fig. 1 • 



The results of our calculation are shown in fig. 2 and 3 for 1.55 21 cm 

wavelengths at nadir observation. For 1,55 cm wavelength, tlm experimental 
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results clearly reflect the surface temperature difference observed durii^ March 
and February flights because at this wavelength the radiation originates from 
very thin surface layer (less than 1 or 2 cm). For 21 cm wavelength on the other 
hand the radiation originates from deeper in the soil and as a result there is 
overlap of the data observed in February and in March. Let us now c<msider 
the effect of surface roughness on the brightness temperature at these two 
wavelengths. On purely physical basis one would expect that as the wavelength 

increases, the effect of surface rcmghness will decrease because at the longer 
wavelengths less of the structure of the reflecting surface will be observed by 
the radiation. The results of our calculation clearly demonstrate this fact. The 
magnitude of the roughness parameter which indicates the structure Observed by 
the radiation is significantly larger for 1.55 cm (h 2 s o.05) than that for 21 cm 
wavelength (h ~ 0.01). The large effect of the surface roughness is easily seen 
in these figures. The flattening of the brightness temperature curve for large 
values of soil moisture cAiserved experimentally can be accounted for the 
surface roughness. 

In fig. 4 we show the data obtained in 1975 and the calculated brightness tem- 
peratures for 21 cm wavelength at nadir observation. The calculation was per- 
formed using the actual ground truth measurement for each test site for which 
the brightness temperatures were observed. Although both observed and calculated 
values are scattered we see that the inclusion of surface rmighness provides good 
qualitative agreement. One should also note in this figure that a linear regression 
analysis of the brightness temperature and the moisture will prodice an excellent 



correlation but the al<^ of the Him d^ncte tq)on the rouf^iness. This dhservation 
has the implication that for any wavelength a lin^r regression analysis may pro- 
dice excellent correlation coefficient but the regression coefficient of that analy- 
sis is not universal in the sense that Uiat it will ^pend up<m the roughness of the 
surface. 

The :»bove discussions clearly demonstrate that for natural agricultural 
terrain the effect of surface roughness on the brightaess temperabire is quite 
significant if die soil moisture is large. Our calculati<Hi indicates that Uiere is 
listening of the brightness temperature for the larger moisture values at all 
wavelengths. Tl» exact moisture value for which the flattening begins depends 
upon the dielectric prq[}erties of the soil at that wavelength. We realize that at 
present we cannot calculate tte value of the r<xighness parameter, h, from a 
first principle consideration and it has been treated as an adjustable parameter. 

Furthermore, the values of the parameter h which give agreement with the 
observed br^htness temperature at 1. 55 and 21 cm wavelengths are not related 
with each other as inverse square of the wavelength which is expected on the 
basis of eqn <9). The surface height variance a calculated using eqn (9) shows 
that the roughness scale important for 21 cm wavelength is in the millimeter 
range and for 1. 55 cm it is of order tenth of millimeter. This variance is the 
same order of magnitude as was observed by Wu and Fung (1972). So that rough- 
ness with amplitudes up to about 1 cm may be the primary cause of the increased 
brightness temperature, and not the corrugations of tte furrowed fields common 
to this area. We are currently formulating the problem so as to model the 
agricultural terrain as a composite rough surface by following the procedure 
given in Wu and Fung (1972). Through such a formulation one may remove the 
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shortcomings this model and provide a better answer to the question tiie 
effect on surface roughness. 

Conclusion 

The problem of microwave emission from a rough half space medium has 
been studied using the radiative transfer equation. 

The dependence of the brightness temperature <m tl» soil moisture of n^ral 
agricultural fields can be explained if one t^es ii^o consideration the surface 
rou^ess. The effect of roughness on the brightness temperature has been found 
to be most severe when the moisture is large. I^cifically it has t^en shown 
that the brightness temperature curve has a natural tendency of flattening when 
the moisture increases to certain value (found to be about 100 percent of field 
caimcity at tiie 1. 55 cm wavelength). The effect of this flattening is to reduce 
the dynamic range of the brightness tem{^rature variation witii soil moisture. 
Further observations of the brightness temperature for h^h values the soil 
moisture is needed to verify the usefulness and the accuracy of the model. 
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lllustrttioiis 

Figure 1. The diel^tric constants for 21 cm and 1.55 cm wavelengths. 

Figure 2. The brightness temperabires for 1.55 cm wavelength at nadir <4)serva- 
ti(m. 

Figure 3. Tte brig^tneis temperatures for 21 cm v^velength at n»iir <^serva- 
tion. 1972 and 1973 data. 

Figure 4. The Bri^tness temperature for 21 cm wavetength at nadir 
observation, 1975 dafo. 
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Figure 2. The brightness temperatures for 1.55 cm wavelength 

at nadir observation. 







Figure 3. The brightness temperatures for 21 cm wavelength 
at nadir observation, 1972 and 1973 data. 
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Figure 4. The brightness temperature for 21 cm wavelength 
at nadir observation, 1975 data. 






